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Abstract 

The earlier introduced method of calculation of quark distributions in hadrons, 
based on QCD sum rules, is improved. The imaginary part of the virtual photon 
forward scattering amplitude on some hadronic current is considered in the case, 
when initial and final virtualities of the current pf, and P2 are different, pf 7^ p^- "^^^ 
operator product expansion (OPE) in pf, P2 is performed. The sum rule for quark 
distribution is obtained using double dispersion representation of the amplitude on 
one side in terms of calculated in QCD OPE and on the other side in terms of 
physical states contributions. Double Borel transformation in pi, is applied to 
the sum rule, killing background non-diagonal transition terms, which deteriorated 
the accuracy in previous calculations. The case of valence quark distribution in 
pion is considered, which was impossible to treat by the previous method. OPE up 
to dimension 6 operators is performed and leading order perturbative corrections 
are accounted. Valence u-quark distribution in vr^ was found at intermediate x, 
0.15 < X < 0.7 and normalization point = 2 GeV^. These results may be used 
as input for evolution equations. 

1 Introduction 



The QCD sum rule approach, invented by Shifman, Vainstein and Zakharov in 1979 [1] 
is now well known as a powerful method, which make possible to calculate in QCD in 
non-model way and with a good accuracy various hadron characteristic like masses, decay 
widths, formfactors etc. The method is based on the operator product expansion (OPE), 
extended to the nonperturbative region. These results were obtained from consideration 
of 2 and 3-point correlators (for review see [0]). A bit later the structure functions - quark 
distribution in photon and hadrons were investigated in the QCD sum rule framework. 
The second moment of photon structure function was considered in |p, and for pion and 
nucleon - in 0, |^, but unfortunately it was difficult to extend this approach for calculating 
higher moments. The general method how to calculate hadron structure functions in 
the region of intermediate x was suggested in and developed in 0. The method is 
based on the consideration of 4-point correlator, corresponding to forward scattering of 
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two currents, one of which has quantum number of hadron of interest, and the other is 
electromagnetic (or weak). In the first order of OPE, in the case, when the hadron is a 
meson, this corresponds to box diagrams hke shown in Fig.l, where p is momentum of 
hadron current and q is momentum of photon. The problem of such diagrams is that 
even if p"^, are large and negative, in the case of forward scattering the singularity 
in t-channel for massless quarks is at t = 0, i.e. large distances in t-channel are of 
importance. However, as was shown in P, 0] the situation changes drastically when the 
imaginary part of the scattering amplitude - the object of interest in case of structure 
functions - is considered. The imaginary part in s-channel (s = (g + pY) of the forward 
scattering amplitude is dominated by small distances contributions at large (negative) 
p^ and intermediate x. (Here the standard notation is used: x is the Bjorken scaling 
variable, x = — g^/2z/, u = pq). The proof of this statement, given in [0], is based on 
the fact, that for the imaginary part of the forward amplitude the position of closest to 
zero singularity in momentum transfer is determined by the boundary of the Mandelstam 
spectral function and given by the equation 

t = -A-^p' (1) 
1 — x 

(it is assumed, that | |^| |)- Therefore, even at t = 0, but not small x and large the 
virtualities of intermediate states in the t-channel are large enough for OPE be available. 
The further procedure is common for QCD sum rule (with some special nuances we will 
discuss later), i.e. dispersion representation on is saturated by physical states and the 
contribution of the lowest particle state is extracted using Borel transformation. In |7]] the 
structure function of nucleon was calculated. Somewhat later, structure function of photon 
has also been calculated p|. But one should note, that sum rule for ci-quark distribution in 
proton obtained in J^, is applicable within rather narrow range of x (0.2 < a; < 0.45) and 
agreement with experiment is not good enough. Moreover, it was found to be impossible 
to calculate structure functions of vr- and p- mesons in this way (that's why the authors 
of P was forced to use special trick, based on VDM, to calculate p-meson structure 
function). The reason for this is that the sum rules, in form used in |^, have a serious 
drawback. 

To understand, what kind of problem it is, let us shortly review the main points 
of the method. Let us consider 4-point correlator with two electromagnetic currents 
and two currents with quantum number of some hadron (for clarity the axial current, 
corresponding to charged pions, will be considered but conclusion is independent of the 
choice of current): 

n^.Ap(pi,P2;gi,g2) = - / e^^^^+^''^™^dWydS(o \T {j,,{x)f;^{y)jr{o)jU^)} I 0) 

jsA = M757AC? (2) 

By considering of forward scattering amplitude in accord with 0, put pi = p2 at the 
very beginning. Among various tensor structures of H/^uXp it is convenient to consider the 
structure (p^iPuPxPp/ ^) ■n(p^, g^, x), and the imaginary part Im H(p^, g^, x) in s-channel is 
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related to pion structure function -F27r(a;)[|. Let us write dispersion relation representation 
of Im n(p^,g^,x) in the variable. As was shown in [0] (see also the correct 

form of dispersion representation is double dispersion relation 

Im n = a{x) + / ^(^l^ + 7 J dudu'p{u,u',x) 



u — p^) J J [u — p'^){u' — p^) 







(We consider lowest twist contributions, the terms of order p^ jq^ are neglected.) In order 
to derive (^) it is convenient to consider first the case, when p\ ^ p\ and go to the limit 
p\^ p\= p^. Then the form (^) is evident. The last term in the right hand side (rhs) of 
(H) represents the propertly double dispersion contribution, the second may be considered 
as subtraction term in variables p\ or p\ and the first term arises as subtraction from the 
second. The interesting for us contribution arises from the pion poles in both variables - 
u and u' in the last term in (^. This term corresponds to the diagram of Fig. 2, where 
the axial current creates the pion, then the process of deep inelastic scattering of virtual 
photon on pion proceeds and finally pion is absorbed by axial current. Evidently this 
term is proportional to the pion structure function. All others in (^) may be considered 
as background. Accept a model of hadronic spectrum, in which p, can be represented 
by contribution of resonance (vr-meson) and continuum (sq is continuum threshold) 

p(m, m', x) = f{x)5{u — m1)5{u' — ml) + p^{x)9{u — sq)6{u' — sq) 

ip{x, u) = (^i{x)5{u - ml) + ip2{x)0{u - Sq) (4) 
where f{x) is proportional to resonance (vr-meson) structure function of interest, 

/(x) ~ 27rF2(x) (5) 

and ipi^2 are some unknown functions, corresponding to non-diagonal transitions. 
The substitution of (^, into (^ gives 

Imtl= ^^^^ +a(x)+ 77 P'i^^^^^')dndu' ^ J ^2{x,u) ^ v^{x,ml) 
(p2 — m^)^ J J {u — p^){u' — p^) J {u — p^) {p^ — m^) 

So ^0 ^0 

The last term in (|D corresponds to Fig.3, where axial current creates a pion, deep inelastic 
scattering proceeds, but the final state is not a pion like in Fig.2, but some excited state 
with pion quantum numbers, which is absorbed by axial current. In order to separate the 
term proportional to the pion structure function - the first term in the rhs of (^, the Borel 
transformation in p^ is applied to (^, which suppresses continuum contributions to (^). 
(The Borel parameter is chosen such that e~*o/*^ ^ 1). After Borel transformation 
we get: 

Bm^ Im ti{p\x) = f{x)^ e-</^'' - y,,{x) e^^H^^ 



^As was mentioned in 0|, the results are more reliable, if invariant amplitude at kinematical structure 
with maximal dimension is used. 
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+ J du J du'p{x,u,u')e-^''+'''^/^^' + J ip2ix,u)e-^ (7) 

So So So 

For the last two terms in rhs of (0), we can assume that p{x, u, u') and ip2{,x, u) are given 
by contribution of bare loop - Fig.l. Because of Borel suppression ~ exp^—So/M"^) <^ 
1 these terms are small and such an approximation does not introduce an essenthial 
error in the final result. However, the second term in rhs of (J^) is not exponentially 
suppressed in comparison with the first. The only way to kill it is to differentiate both 
sides of (multiplied by exp^ml/M^)) over Just this procedure was used in 

^ to determine nucleon structure functions. But, as is well known, the differention of 
approximate relation may seriously deteriorate the accuracy of the results. In QCD sum 
rules such procedure increases contribution of nonperturbative corrections and continuum 
contributions, sum rules become much worse or even fails (as for p-meson). For vr-meson 
situation is even worse, because direct calculations show, that bare loop contribution 
corresponds only to non-diagonal transitions. 

In this work we suggest the modified method of calculation of the hadron structure 
function, which is free from this problems and is completely based on QCD sum rules. 
We will illustrate it on an example on the vr-meson structure function calculation, which 
usually is much "dangerous" case. 



2 The idea of the method. 

The idea of the method is to consider at the begining non-equal pi ^ p\ in (|^) and perform 
all calculations for this case. Instead (|^) dispersion representation takes the form 



V9(x,m) 



u 



Pi 



du + 



!f{x,u) 



u 



du+ dui / du 



p{x,Ui,U2) 



'{Ui -Pi)(m2 - pi) 



Apply to d^) double Borel transformation in p1,P2- This transformation kills three first 
terms in rhs of (|^) and we have 



^M^^M^Im U{pI,pI,x) 



dui / du2p{x,Ui,U2)exp 



Ui 



Ml 



U2 

Ml 



(9) 



One can divide the integration region over -Ui 2 into 4 areas (Fig.4): 

I Ml < So;U2 < So] 

II Ui < So; U2 > So 

III U2 < So; ui > So 

IV Mi,2 > So 

Using the standard QCD sum rule model of hadronic spectrum and the hypothesis of 
quark-hadronic duality, i.e. the model with one lower resonance plus continuum, one can 
easily notice]^, that area (I) corresponds to resonance region. Spectral density can be 
written in this area as 



^We restrict ourselves by simplest model, because higher resonance contribution in any case will be 
suppressed after double borelization. 
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M2, x) = fl- 2itF2{x)6{ui - ml)5{u2 - ml) (10) 

where is defined as 

(0 I JA5 I tt) = Upx U = 131 MeV 

In area (IV), where both variables ui^2 are far from resonance region, the non-perturbative 
effects may be neglected, and, as usual in sum rules, spectral function of hadron state is 
described by the bare loop spectral function in the same region 

P(Mi, U2, x) = p°(Ui, U2, X) (11) 

In the areas (II), (III) one of variables is far from resonance region, but other is in the 
resonance region, and spectral function in this region is some unknown function p = 
'ip{ui,U2,x), which corresponds to transitions like it continuum, as shown in Fig. 3. 
After double Borel transformation the total answer for physical part can be written as 
{Ml, M| are Borel masses square) 



BiB2[Im U] = 27rF2{x) ■ f^e ^ + dm du2ip{ui,U2,x)e ^ ^ 

so 

+ [ dui [ du2ij{uuU2,x)e~^'^^^\ [ [ duidu2p\ui,U2,x)e~^^^'^'' (12) 



so So So 

In what follows we put for simplicity = M| = 2M^. The one of advantages of this 
method is that after double Borel transformation unknown contribution of (II), (III) areas 
(second and third term in (|l^)) are exponentially suppressed. Using duality arguments we 
estimate the contribution of all non-resonance region (i.e areas II, III, IV) as contribution 
of bare loop in the same region and demand their value to be small (less than 30%). So, 
equating physical and QCD representation of 11, and taking in account cancellation of 
appropriate parts in left and right sides one can write the following sum rules (we omit 
all terms, which are suppressed after Borel transformation) 

Im ^QCD + Power correction = 27rF2(x)/^ 

SQ So 

^mIlQCD = J J p {ui,U2,x)e ^S^; (13) 



(The pion mass is neglected.) It can be shown (see Appendix), that for box diagram 
ip{ui,U2, x) ~ S{ui — U2) and, as a consequence, the second and third terms in (P^) are 
zero in our model of hadronic spectrum. 
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3 Calculation of box diagram. 



The diagrams, corresponding to unit operator contribution, are shown in Fig.la,b. Note, 
that crossing diagram, shown in Fig.lc does not contribute, their contribution found to 
be in leading twist. (This is a sequence of kinematics, so such crossing diagrams also 
are zero for higher dimension corrections in the leading twist.) 

It is enough for us to calculate the distribution of valence w-quarks in pion, since 
d{x) = u{x). For this reason restrict ourselves by calculation of Im 11 for the diagram 
Fig. la. 

Consider first the case pi = p2 and demonstrate, as was announced in Sec. 2, that in 
this case the contribution of box diagram attributes only to non-diagonal transitions, like 
in Fig. 3 and refers to background terms in (|^). Diagram Fig. la contribution is equal 

/mn,.,. = / {k + qm {p - k-Y] X 

X Tr[ -fxhX^^ + q)iyha{k - p) ] (14) 

Calculate the trace and omit the terms, which cannot contribute to the interesting for us 
structure ~ PfiPuPxPa/^- We get 

12 r d^k 

Im n^,,. = --2] jr^kukxik. - PaMk + qYmp - kf] (is) 

Calculation of the integral leads to: 

3 1 

Im U^iyXcT = P^PuPxPa — ^a;^(l - x) (16) 

(only the terms ~ p^PuPxPa are kept) and 

Jm n(p^x) = -^^x^(l -x) (17) 
Substitute ( |T7| ) into (^ and perform Borel transformation. We get: 

l^^il _ _ e-^o/M^) = 27r/>.(x)^ + Mx), (18) 

where Ut^Ix) is the distribution of valence u quarks in pion (pion mass is neglected). 
Looking at dependence in (|18D it becomes evident, that in this appropach the attempt 
to separate the pion contribution from the background by studying dependence (e.g. 
differentiation over 1/M^) is useless - up to small correction ~ g-^o/M ^j^^ ^^^^ diagram 
contributes to the background only. 

Consider now the more promisable approach, pf p\- Since nonequality of p\^ p\ 
is important for us only for Borel transformation, i.e. in the denominators of dispersion 
representation (^), in the calculation of numerator, resulting in kinematical structure 
PiiPvPxPa we can put p\= P2= P- Therefore, in order to understand the essential features 
of corresponding integrals in case of non-equal pfyp"^, it is sufficient to study insread of 
(|T^ a more simple integral 
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ImT{plplq\u)= fd'k^^-—^ S[{k + qmip,~kr] (19) 

The direct calculation of the integral in rhs of ([ToD (see Appendix) shows, that it may be 
represented in the form 

2v/x 

ImT{pl,pl,q'^,u) = - — / 2 2^" (20) 

Aux J u~ pfu — P2 

(Higher order terms in pl/q"^, pi/^f are neglected.) At p\ = p^ it gives 

Im np\ q\ y) = ^2 (21) 

as it should be. (pO]) may be rewritten in the form of double dispersion representation (H) 
with a{x) = fi^x) = and p{u,u',x) proportional to S{u — u') 

2 2 \ ''^ f ^{"^ — u') , , / 
ulm T{p,,p2,x) = — — / 7 TTrdu du (22) 

(Higher twist terms are omitted). From this consideration it becomes clear, that in order 
to go from the case of p\= p\ = p^ in the calculation of the box diagram Fig. la (|14D to 
Pi 7^ pi 5 it is enough to substitute in the final result the factor Xj-p? by0 

Therefore instead of (0) we get 

oo oo [■/ /\ 

A ^ -Al-.)liul,u'^-^Z^ (24) 

Perform double Borel transformation in p\.,p\. It kills nondesirable depending on one 
variable subtraction terms in (^ and we have the sum rule for valence w-quark distribution 
in pion 

where it was put Mf = M| = 2M^. (As is known [11] the charactiristic values of Borel 
parameters Mf , M| in double Borel transformation are about twice of Borel parameters 
in ordinary Borel transformation, used in mass calculations.) 

Before going to more accurate consideration with account of higher dimension op- 
erators and leading order (LO) perturbative corrections, let discuss the unit operator 
contribution in order to estimate, if it is reasonable. The calculation of the pion decay 
constant f.^, performed in in the same approximation results in 

■^It must be mentioned, that such substitution is vahd only for box diagram, it does not take place for 
more complicated diagrams, considered in next Section. 



7 



/^2 = _i_M^(l-e--/^^^) (26) 



Substitution ( p6[ ) into (|2^) gives 



u^{x) = 6x{l - x) (27) 

One can note, that 

1 

J u^{x)dx = 1 (28) 



in agreement with the fact, that in the quark-parton model there is one valence quark in 
pion. Also, it can be is easily verified, that 

1 

J xu^{x)dx = 1/2 (29) 



what corresponds to naive quark model, where no sea quarks exist. So one can say, that 
formally unit operator contribution corresponds to naive parton model. 

Of course, eq.(P5|) has only formal sense, because, as was discussed in Introduction, 
our approach is correct only in some intermediate region of x. The boundaries of x, 
where this approach is correct, will be found, if one takes into account nonperturbative 
power corrections. In the next section we will discuss them. At the end of this section let 
us discuss perturbative corrections. We take into account only LO terms, proportional 
to Zn(Q^/M^), and choose = — 2GeV^ - for the point we calculate our sum 
rules. Finally, the result for bare loop has the form (the second term in square brackets 
corresponds to perturbative correction, taken into account). 



u^[x 



3M^x 
271^ P * 



-x] 



14 



[l/x+4:ln{l-x)- 



2(1 - 2x)lnx^ 



X 



'1-e 



-So/M2 



In the calculation we choose the normalization point = /i^ 



(30) 

The fact that we take 



into accounrt as-correction at the point = 2 GeV^ means that our final results for the 
structure function (we write it in the next section) can be used as an input for evolution, 
starting from this value of Qq. 



4 Calculations of higher order terms in OPE. 

In this section we discuss the power correction contribution to sum rules. The power 
correction with lower dimension is proportional to gluon condensate (G^^G^^) with d = A. 
As was discussed above, only s-channel diagrams (Fig. la) exists in the case of double 
borelization. (G^^G^^) correction was calculated in a standard way in the Fock-Schwinger 
gauge x^A^ = 
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The quark propagator iS{x,y) = {ip{x)ip{y)) in the external field has the well 
known form |jl3[ (our sign of g is opposite to that of |[T3|1 ): 



iS{x, y) = iS^{x — y) — g J d^ziS^{x — z) ■ iA{z)iS^{z — y) + 

g^ J J d^zd^z' iS\x - z)iA{z)iS\z - z') ■ %A{^) ■ zS°(^' - 2) + ... (31) 
Here 5*° is free quark propagator; A = (l/2)A"7^y4'^ and 



A^ix) = i-)x,G;^ + {-)x^x,[D^G;^m + i-)x^x,xp[D^DpG;^] (32) 



When calculating one should take into account quark propagator expansion up to the 
third term and only the first term in the expansion of the external field (Fig. 5). 

These diagrams have been calculated using the program of analytical calculation RE- 
DUCE. Surprisingly, in the case of the double borelization the sum of all diagrams Fig. 5 
was found to be 0. So, the gluon condensate contribution to the sum rule is absent. 

Before we discuss the d = 6 contribution, let us make the following remark. Due 
to the fact that we are interested only in the intermediate values of x, we should take 
into account only loop diagrams. Really, one can easily see, that the diagrams with no 
loops (like those in Fig.6) are proportional to 6{1 — x) and is out of the region of the 
method applicability. There are a large number of loop diagrams, corresponding to = 6 
corrections. First of all, there are diagrams which correspond to interaction only with 
gluon vacuum field, i.e. only with external soft gluon lines (see Fig. 7). Such diagrams 
may appear, if we take: 



a) all possible combinations, which appear when expansion of quark propagator (pl|) 
is taken into account up to the fourth term and in expansion of the external gluon field 
( ^21) only the first term is kept. For example, it is the fourth term of the expansion for one 
quark propagator and the first term (free propagator) for other three (Fig. 7a), the second 
term of the expansion for three quark propagator and one propagator is free (Fig. 7b), 
the third term of the expansion for one quark propagator and the second term for other 
(Fig. 7c) etc. 

b) all possible combination, when the second and the third terms of expansion of gluon 
field (|3^) is taken into account, like those, shown in Fig.8. 



The diagrams of Fig. 7 are, obviously, proportional to {g^ f°'^'^G'tyG^^(^G'p^) and when 



^ fiu^ a(5^ pa 

calculating it is convenient to use representation of this tensor structure given in 



(0 I g' r'^Gl^Gi^G% I 0) = (1/24) (0 | r'^Gt^filGl^ \ 0) ■ {g,Mpp+ 
+g^ii3gapgav + gaag^ipgup + gpvgp.agi3a— 
~ gppgaagpu — gp^agu^gap — gaugppg^a — g^pg^ag va ) (33) 

The diagrams of Fig.8 are proportional to (0 | ^pG^'^^D^-G"^^ \ 0) and 
(0 I G^-^^D pD tG°^^ I 0) . Using the equation of motion it was found in |^ 
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20" 



-(0 I D,G%D,Gl^ I 0) = (0 I D,G%D,D^Glp \ 0) 



p) + 2^9fil39aa9pu + 9au9pp9l3a " 9aa9pp9ul3 " 9pu9pa9l3a) 



+ 



O^ = -{0 I qXj; I 0) ± ^(0 I 9r''G%Gl,Gl^ \ 0) (34) 



a9au9l3p + 9p,l39ap9cTU 9pcT9ap9ul3 9pl39na9v 

where = J2 tpilp^ ^i- 

i 

From (p3|), ( |34D one may note, that these tensor structure are proportional to two 
different vacuum averages: 



fjl\0) and {0 \ g-'G^^Gl^G^J^'^ 



First of them (0 | g'^j'^ \ 0) by use of the factorization hypothesis easily reduces to (giptp)"^, 
which is well known. 

(0|/j;|0) = -(4/3)[(0|#^|0)]2 (35) 

But (0 I g'^G'I^^GlpGp^f"'^'^ I 0) is not well known, there are only some estimations based 
on the instanton model |T^. Fortunately, in the sum of all diagrams of this two types 
all terms proportional to this dimension 6 gluonic condensate are exactly cancelled, and 
the sum of diagrams of Fig. 7 and Fig. 8 is proportional only to {gipip)"^. 

We consider now an another type of diagrams which also give contribution to = 6 
power corrections. Such diagrams appear when the external quark field is taken into 
account, i.e., one should take into consideration the expansion of quark field: 

V'(x) = z^(0) + x^, [V„,V^(0)] + ]<t^,x^, [Va, V,,V^(0)] + ... (36) 

where V is co variant derivative. 

In this case there appear diagrams like those in Figs. 9-11, where quark (and antiquark) 
line is expanded and the first and the second terms of the expansion ( pBf ) are taken. The 
expansion of the external gluon field (^) is also accounted up to the second term. For 
the diagrams of Fig. 10 gluon propagator in the external field is also accounted (we discuss 
it a bit later). 

All these diagrams can be divided into two types with quite a different physical sense. 
The first type of diagrams - like those in Fig. 9 - corresponds to the case, where all 
interactions with vacuum proceeds out of the loop. Such diagrams correspond to logarithm 
corrections (evolution) to the corresponding non-loop diagrams (without hard gluon line). 
Since, as was discussed in Sec. 3, we will not take into account these non-loop diagrams, 
then it seems reasonable, that at the same level of accuracy we do not take into account 
their evolution. So, all the diagrams of this type should be omitted. The problem of 
correct calculation of non-loop diagrams and their leading logarithmic correction is a 
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special problem, which will not be discussed here. In any case, estimations and physical 
reason show, that their contribution would be significant at large x and negligible in 
intermediate region. We shall see at the end of the paper, that sum rules themeselves 
indicate region of x where effects of the non-loop diagrams and their evolution may be 
neglected. 

So, according to the previous discussion, we should bear in mind only those diagrams, 
where interaction with vacuum takes place inside the loop. (Figs. 10, 11). Such diagrams 
cannot be treated as evolution of any non-loop diagrams and are pure power correction 
of dimension 6. All these diagrams are, obviously, proportional to 

(0 I ^'pDfil, I 0); (0 I ^^(V.V^J)G;:, I 0); 



(0 I {VM)i^'pGl, I 0) 

These tensor structures were considered in where using the equation of motion the 
following results have been obtained 



\ / /3o 



n\bd 



(0 I (Vai'aY'ip^pG'^j^ I 0) can be easily calculated using the results of p 
For diagrams in Fig. 10 we use the following expansion of gluon propagator 



—7 r fl'^k 

^ rnpl / "* -iku 
J 



-ikxVaG^x - T^Kyay/skx- 



'^{kH^x - 2k^kx)){D^Gpxy + l(yc^ + ^) {DxG^xY 



5up + 2 



G''^p + 2i—{DaGupf' 
(38) 

This expression can be found using the a method of calculation of gluon propagator 



in external vacuum gluon field, suggested in [|12[. The same result up to ~ G term is 
explicitly written in [1^ (see also Jl^). The total number of c? = 6 diagrams is enormous 
— about 500. All of them were calculated by use of REDUCE program. The final result 
for d=6 correction after double Borel transformation have the form 



Im 11'^=^ 



1 2 {ga? 1 



(27r) 



rg ■- 



M4 36-25 



+ 4740x^ + 8847x^ + 2066x^ - 2553x + 1416 



-5784x^-1140x^-20196x2+20628x-8292) ln{2) 

1 



xil — x) 



(39) 
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{ga)"^ = Anas ■ (27r)'^(0 | ipip \ 0)^. Before we write the final result of the sum rules, let us 
make one note. One can see, that in contribution of d = 6 operator (p9D strong coupling 
constant appears as factor, and again it appears in structures {gaY- The factor g"^ 
corresponds to interaction with quark propagator (vertices of hard gluon line in diagrams 
in Fig.9,10, or vertices of external gluon in diagrams in Fig.6,7), and it is reasonable to 
take it at the renormalization point /i^ = Q^. On the other side, g'^ in structure {gaY 
appears as a consequence of use of equation of motion and its normalization point should 
be taken in such a way that the quantity ^^(O | ipip \ 0)^ is renormalization group invariant. 
Finally, substituting results for bare loop (0) and power corrections (|39|) , we can write 
the sum rule for quark distribution function in pion: 



XUt,{x) 



3 A/f^ 



x^{l-x)- 



1 + 



'a,(M2) • ln{Ql/A'P) \ ( I + Ax ln{l - x) 2(1 - 2x)ln x" 



37r 



X 



X 



) 



X] 



(40) 



(27r)4 . 37 • 26 . M6 x3(l 

where uj{x) is the expression in square brackets in (^). We choose the effective LO QCD 
parameter Aqcd = 200 MeV, Qq = 2 GeV^. The value of renorminvariant parameter is 
equal 



asa^ = a,(M2 = 1 GeV'^) ■ (0.55 GeV^ = 0.13 GeV^ 



The value of a was taken from the best fit of the sum rule of nucleon masses (see 
Appendix B). Continuum threshold was varied in the interval, 0.8 < Sq < 1.2 GeV^ and 
it was found that the results only slightly depend on it. The analysis of the sum rule (^) 
shows, that they are fulfilled in the region 0.15 < x < 0.7; power corrections are less than 
30%, and continuum contribution is small (< 25%). Stability in Borel mass parameter 

dependence in the region 0.4 GeV^ < < 0.6 GeV^ is good; especially in the region 
of X < 0.4 the dependence is almost constant (see Fig. 12). 

The final result for u^{x), (at = 0.45 GeV^, So = 0.8 GeV^), is shown in Fig.l3 
(thick solid line). On Fig. 13 is also plotted the curve of w-quark distribution in pion, 
found in by using evolution equation and the fit to the data on Drell-Yan process, 
performed in pO|). Bearing in mind, that NLO a^-corrections are not accounted and one 



may expect, that they would increase Ut,{x) at low x and decrease at large x, one may 
consider the agreement as good. We also show in the same figure pure bare loop contri- 
bution (line with squares) and contribution ( |5(]| ) of bare loop with perturbative correction 
(crossed line). One can see, that pure bare loop is not in a quite good agreement with 
experiment and both perturbative correction and power correction improve the agreement 
with experiment. Let us discuss, why stability became worse when x became larger (see 
Fig. 12). From our point of view, it reflect the influence of non-loop diagrams (and their 
evolution), which were not accounted as it was discussed in sect. 4. Indeed, the non-loop 
diagram which formally are proportional to 6{1 — x), of course really would correspond to 
some function with maximum close to x=l and fast decreasing when x decreased. That 
is why effects of such diagrams (and their evolution) are negligible at x < 0.4 — 0.5, but 
may be more or less sensible at large x, and deterioration of stability probably reflects 
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this fact. We repeat, finally that obtained valence w-quark distribution function Ut^{x) 
can be used as input for evolution equation (starting from point Ql = 2 GeV"^). 

Let us now discuss at the end the estimations for the moments of quark distribution 
which can be found with the help of the results obtained. 

To get the moments, one should make some suggestions about the region of small 
X < 0.15 and large x > 0.7 0.8, where our method is inapplicable. Of course, in this 
case the estimation of moments are not model-independent and the accuracy of estimation 
of moments should be treated as lower than for the structure function (^) itself. 

If we make a natural supposition, that at x < 0.15 Ut,{x) ~ according to Regge 

behaviour, and at large x > 0.7, u.j,{x) ~ (1 — x)^ according to quark counting rules. 



then, matching these functions with our result (|40|), one may find, that 



1 



Mo = J u.„{x)dx ^ 0.84 



Ml = J xuT,{x)dx ^ 0.21 



at = 0.4-^0.45Ge^^. These results only slightly depend on the choice of the points of 
matching (not more than 5% when we vary lower matching point in the region 0.15 ^ 0.2 
and the upper one in the region 0.65-r-0.75. One may note that Mo which has the physical 
meaning of the number of quarks in pion (and should be Mq = 1) is really close to 1 
within our accuracy ~ 10 20%. Mi has physical meaning of the part of momentum 
carried by a valence quark, and the value Mi ~ 0.21 is in good agreement with well 
known fact, that two valence quarks in pion take about 40% of the total momentum. 
This work was supported in part by RFBR grant 97-02-16131. 
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Appendix 

In this Appendix the double dispersion representation (^) of integral ([T9| ) is proved. 
It is convenient to change variables in ([T9|) and put p — k = k'. Then ([T9|) takes the form 
(prime is omitted) 

ImTiplplq',u)= ld'kj-^j-^5[ip + q-kY]S{k') (Al) 

Assume that = q'"^-, t = (pi — ^2)^ = and choose the Lorenz system, where 4- vector 
P = [Pi +P2)/2 has only 2;-component equal P^. From 

t=iPi- P2f =pi+pl + 2piP2 = (A.2) 

it follows 

P' = -P^ = {pI+pI)/2 (A.3) 
Introduce 4-null- vector r = pi — p2, r'^ = 



We have 



Use the notation 



Then 



rP=ipi-pi)/2 (A.4) 



p^ = P + r/2, p2 = P-r/2 (A.5) 



qP = v = qpi = qp2 (A. 6) 



q'pi = qp2 = I' + {p\- pi)/ 2-, qr = q'r = (A. 7) 



We can choose the coordinate system, where 4-vector q has only time and 2;-components 
and 



= -i//v^=P2, qo = Jv^ - q^PyV^ ^ u/V^ = -q, (A.t 



The last equality corresponds to account of lower twist terms. From ( |A.4| ) and ( |A.7| ) for 
4-vector r with components r = {tq, r±, r^} it follows 



1 pj-pl 1 1 pi-pI I pI-pI 

Tn = —U — — ^^^^^^^^= ^ — 7^ = — — ~ 

2 sT-P^ sfU^^^q^ 2 ' 2 °' 

= _i ,^M:i|£ (A.9) 



The components tq and r^. are equal in the lowest twist approximation and of order y/—P'^ 
if pI ~ P2, while r± ~ (p^ —p'i)!^^^'^^ i-^- of next order in this approximation and may 
be neglected. The argument of the first 5-function in ( [A.l| ) is equal 
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ko + r±k±Cosip = 



-p2 ' 2 

(A.IO) 

where is the azimutal angle between p± and k±. The the last term in (|A.1CI|) may be 
omitted - it is of the next order in p^/ z/: it may appear only squared because of integration 
over (fi in (|A.1|) ). (This fact can be proved by direct calculation.) From the inequality 



kl - kl > 0, 



(A.ll) 

the inequality follows, which defines the integration domain over k^ in the integral ([A.l|) : 



1 



It is convenient to use the notation 



The integration domain is 



V 



2V^k, 



(A.12) 



(A.13) 



-2u <v < -P^{l-x) 
The denominators in (|A.1| ) are calculated by using the relations: 



(A.14) 



(p^ - kf =p{- 2pik = pI- 2Pk - rk (pa - k^ = pj - 2Pk + rk 



(A.15) 



Pk = -y^^k,, rk^ro{ko + k,) =roV^{l-x) = -{pI-pI){1-x) (A.16) 



(In the above equalities ( |A.1(]| ) and ( |A.9D were exploited). As a result we get ((5-functions 



were eliminated by integration over k'j_ and ko): 

-P\l-x) 

TT f 1 

Im T = , / dv 



4goy^ 4 ^''pl + v-ipl-pl)il~x)/2 pl + v + ipl-p^2)i^-x)/2 

(A.17) 



Changing variables 

t; = -p2(l -X) -MX (A.18) 

gives the final answer 

TT 7 1 1 

ImT = / ^ ^du (A.19) 

Avx J u — pfu — P2 

(The upper limit of integration was put as infinity, what is legitimate in the lowest twist 
approximation) . 
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Figure 1: Diagrams, corresponding to unit operator contribution. Dashed lines with 
arrows correspond to photon, thick sohd - to hadron current 
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Figure 2: Diagram of forward photon-pion scattering. 
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Figure 3: Example of non-diagonal transition 
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Figure 4: Integration region in double dispersion representation. 
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+ all combinations 



Figure 5: Diagrams, corresponding to d = 4 operator contribution. Dashed lines with 
aroows correspond to photon, thick sohd - hadron current, wave hues correspond to 
external gluon field 
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J V 




Figure 6: Examples of non-loop diagrams of dimension 4. Wave lines corespond to gluons, 
dot means derivative, other notations as in Fig.l 
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Figure 7: Diagrams of dimension 6, see text. All notations as in Fig.6 
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+ all combinations 



Figure 8: Diagrams of dimension 6. External gluon line with dot corresponds to deriva- 
tives of gluon lines. All notations as in Fig.6 
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Figure 9: Dimension 6 giagrams with out the loop vacuum interaction. All notations as 
in Fig.8 
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Jl combinations 



Figure 10: Diagrams of dimension 6, corresponding to quark propogator expansion 
(eq.37). All notations as in Fig. 8 
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all conibiriaLions 



Figure 11: Diagrams of dimension 6, corresponding to quark and gluon propogator ex- 
pansion (eq.37, 38). All notations as in Fig.8 
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Figure 12: Borel mass dependence of quark distribution function at various x 
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xu(x) 




Figure 13: Quark distribution function in pion, noted "total". For comparison fit from 
[19], noted "GR", is shown. Also bare loop ("bare") and bare loop with perturbative 
corrections (noted "1"), are shown 
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